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ABSTRACT 

Weak lensing alters the size of images with a similar magnitude to the distortion due 
to shear. Galaxy size probes the convergence field, and shape the shear field, both 
of which contain cosmological information. We show the gains expected in the Dark 
Energy Figure of Merit if galaxy size information is used in combination with galaxy 
shape. In any normal analysis of cosmic shear, galaxy sizes are also studied, so this 
is extra statistical information comes for free and is currently unused. There are two 
main results in this letter: firstly, we show that size measurement can be made uncor- 
related with ellipticity measurement, thus allowing the full statistical gain from the 
combination, provided that -s/Area is used as a size indicator; secondly, as a proof of 
concept, we show that when the relevant modes are noise-dominated, as is the norm 
for lensing surveys, the gains are substantial, with improvements of about 68% in the 
Figure of Merit expected when systematic errors are ignored. An approximate treat- 
ment of such systematics such as intrinsic alignments and size-magnitude correlations 
respectively suggests that a much better improvement in the Dark Energy Figure of 
Merit of even a factor of ~ 4 may be achieved. 
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1 INTRODUCTION 

Weak gravitational lensing by the intervening nonuniform 
matter distribution has been recognised as a potentially 
very powerful tool for probing the growth rate of poten- 
tial fluctuations and the geometry of the Universe through 
the distance-redshift relation. Traditionally the statistic of 
choice has been cosmic shear — the distortion in the shape 



of the image of a source (see Munshi et al. 2008 and ref- 
erences therein). However, weak lensing has other effects, 
such as a magnification of the size of the image, and a corre- 
sponding change in the flux of sources. In an ideal analysis, 
one would like to use all of this information. Flux magniflca- 
tion is beginning to be explored ( van Waerbeke|20l0 Hilde- 



brandt et al. 



2009 



2013 Duncan et al. in prep) and after an 
early study ( Bartelmann et aL]|1996[ ) size magniflcation has 
begun to receive attention, both theoretically (Casaponsa 



et al.|20l2 |, and observationally ( Schmidt et al.|2012 l. Cas- 
aponsa et al.| ( |2012[ ) showed that the convergence fleld can 
be recovered from the measured sizes of simulated galaxy 
images without any evidence of bias, provided the galaxies 
are larger than the point-spread function (PSF) and have 
S/N larger than 10. These are very similar requirements for 
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accurate estimation of shear, and since the shape measure- 
ment process also inevitably investigates size, this informa- 
tion comes for free. The focus of this letter is two-fold: flrstly 
to analyse what extra information is provided by size, and 
secondly to demonstrate that size and shape measurements 
can, with a careful definition of the size, be made uncorre- 
lated, so we can use the full statistical power from adding 
size measurements. It is intended to be the second step in a 
programme to develop more powerful combinations of weak 
lensing measurements to extract the full statistical power, 
and a number of questions are not addressed in this study, 
whose purpose is a proof-of-concept to illustrate that signif- 
icant gains are possible. With reasonable assumptions, we 
find that Figures of Merit for Dark Energy studies may be 
improved by significant factors, with no additional observa- 
tional data required. 



2 STATISTICS OF COMBINED SIZE AND 
SHEAR MEASUREMENT 

In this section, we study what improvements in error bars 
we might expect from combining measurements of size and 
shape. As the result is not quite as one might expect, we first 
illustrate the effect with a simplified case (a single source 
plane and single mode) , before performing Fisher matrix cal- 
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culations to analyse the efTect on a future survey designed 
to produce a large Dark Energy figure-of-merit. We ignore 
systematic effects in this section, and consider them in Sec- 
tion H 

Lensing effects are described by the transformation ma- 
trix mapping source angular positions to image positions, 



K - 7i 
-72 



1 



—72 
K + 71 



which defines the convergence field k and complex shear field 
7 = 71 -f 172. The magnification of surface area elements, v 
is given by the determinant: 



(1) 



If |k| and I7I ^ 1 (which we assume throughout) this can 
be approximated by ~ 1 + 2k, so a length scale defined 
by the square root of the area, which we will see is a very 
useful definition of size, will scale to linear order by 1 -I- k. 



In the Limber ( 1954 1 approximation, the angular power 



spectrum of the lensing potential between tomographic red- 



shift bins i and j is given by ( Takada fc Jain(2004J 
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where x is comoving distance, P{k; x) is the 3D matter 
power spectrum and Xm(x) is the transverse comoving dis- 
tance corresponding to comoving distance X- The lensing 
weight functions u;(i)(x) are given by 



W(»)(x) 



r^+i , dx p(x) ^"^^'"^^ 
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(3) 



where p(x)'^X — n(z)dz for a galaxy redshift distribution 
n{z), Ui is the number density in the zth bin, Xi ^^'^ Xi+i 
are the boundaries of the ith tomographic bin. On the full 
2D sky, the spherical harmonic expansion coefficients of the 
shear and convergence fields (associated with a particular 
tomographic bin) are related to those of the lensing potential 
(e.g., , Hu|2000} [Castro et"aL||2005| ) : 



1 1 2 
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Taking as the estimators k — ln(A/A) = k + ln(As/A) and 
7 = e = 7 -I- Bs , where A is the mean size at the appropriate 
redshift and As is the unmagnified source size, the (cross) 



power spectra are given by (e.g. Hu (2002l) 



A7271 
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(5) 



where ae and crinA are the dispersions in the intrinsic (com- 
plex) ellipticity and log-size InAs respectively. 

Since K^m, 7i, «m and 72, are complex, care must 
be taken in constructing the covariance matrix to en- 
sure that all of the information has been included cor- 
rectly. Here we take our data vector to contain entries 
for the expansion coefficients and their complex conju- 

')^ for the combined shear- 



)""" for the shear only 
and 



gates, i.e., d'"'^) = (z'^-t', z(«'T> 
convergence data and d'''' — (z'^' 
case, where z'"'^' = (... 71, 72, ... ) 

z^^-* = (. . . 71, ^m(i), 72, em{i) ■■ ■)'^ contain the fuU set of rel- 
evant complex coefficients. Note that populating the data 
vector with the real and imaginary parts of the relevant 
fields explicitly is entirely equivalent (see, e.g., |Picinbono 



1996 1. Furthermore, to avoid duplication of information only 



m ^ modes are included and care must be taken not to 
double count the m — modes for which neo, 71,^0 and 72,^0 
are real (recall that since is a real field, (jjim oc (jji-m ''■nd 
(jieo G K). The full covariance matrix T of the data is defined 
as: 

(zzt) 
(zz)* 

C 
C 



r = (ddt) = 



(zz) 
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where in the second line we have used the fact that (zz) — 
and C = (zz^) G R. Since different £ and m modes are un- 
correlated for an all-sky survey (we relax this later), C will 
be block diagonal with each {£, m)-mode contributing one 
diagonal block: 



where 



I -I- n 



n = diag(ni, 712 ... ), 
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and (8> is the tensor product. 

The Fisher matrix, Fca, is the negative expectation of 
the second derivative of the log-likelihood with respect to the 
model parameters labelled by a and p. If the data can be 
assumed to be Gaussian distributed with fixed means, such 
that the covariance matrix is determined by the parameters 
of interest, the Fisher matrix can be computed from the 
covariance matrix and its derivatives ( Tegmark et al.|1997 |: 



^Trace [r^'r, ^.F^'r, a] 



Trace fC^^C, ^C^^C 
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where a subscripted comma refers to derivatives with re- 
spect to the foUowing parameter. Since C is block diagonal, 
with each (l,m) mode contributing a block Ci, the Fisher 
matrix can be written as a sum over modes 



t, /3J 



(10) 



where we have also included a factor /sky to approximately 
account for incomplete sky coverage. To illustrate, let us 
consider estimating the amplitude of the lensing potential 
power spectrum Ci from a single mode and a single topo- 
graphic bin (so we drop the (i) subscript), ignoring for now 
systematic effects in shape and size, and assuming the same 
number density for both size and ellipticity. In this case, for 
a shape and size analysis, the covariance matrix C, Eq. ([7| 
is 



i(«;,7) 



Cl 



Ct 
Ct 

Cl + O-lnA/^i 



(11) 



where ct = t^Ct/A. For a shape-only analysis, C is the top- 
left 2x2 sub matrix. The Fisher matrix given by Eq. Q 
is a scalar in these cases, and reduces to (multiplying by 2 
since the covariance matrix has two blocks of C), 



(12) 



and 



in' 



{2nci + a'iy 

Therefore the error bar on ct is reduced by a factor 

{2R+1){2S + 1) 



•^shca 



•^shcar -|- size 

2{2SR + S + R) 



(13) 



where R = crfnA/'^e and S = nci/a'^ is a measure of signal- 
to-noise. 

We see that in the high S/N limit, there is no gain; es- 
sentially both size and shape are measuring the same quan- 
tity with a vanishingly small error bar. Since the signal we 
are using here is the variance around the zero mean, there 
is no benefit. The other limit is interesting; in the low S/N 
regime, the gain is a factor of 1.5 if = o^i^x- i.e. we esti- 
mate the variance with an error smaller by a factor 3/2. 



3 UNCORRELATED AREA AND SHAPE 
MEASUREMENT 

We have so far assumed that the estimates of the shear and 
convergence are uncorrelated. It is not obvious that this 
case be achieved, but in this section we demonstrate that 
for galaxies which have exponential brightness profiles, the 
estimate of y^area is uncorrelated with the estimate of ellip- 
ticity when estimated using model-fitting methods such as 
lens&t ( [Miller et al.|2007[ ). For more complex morphologies, 
we would expect the correlations to be non-zero, but there 
is a reasonable expectation that they would be small. 

We model the galaxy surface brightness fi with a thin, 
intrinsically circular, disk with an exponential profile with 
scale length R. The apparent shape of the galaxy image is 
determined by the angle 77 between the disk normal and the 



line of sight. If the projected elliptical image has a position 
angle (j), then after some algebra, the surface brightness may 
be written, as a function of polar coordinates (r, ^p), 



^(r, Tp\X, e, (f), fio) = Mo exp 



2ecos2(V' - (f>) 



where is the central surface brightness. 



(14) 



A : 



1 + e 



and a, b are the semi-major and semi-minor axes, e is the 
(magnitude of the) ellipticity, defined by e = {a — b) /{a + b). 
The area of the ellipse is nab = ttA^, so A is a measure of 
the square root of the area. We ignore errors in the centroid 
in what follows. 

If we estimate the four parameters of the model, — 
(A, e, ^0, (j)) from a set of pixels with gaussian white noise er- 
rors, then we can obtain the covariance of the estimates from 
the Fisher matrix, here with the covariance matrix fixed 
(e.g., [Tteimark et al.|19"97t : 



:Tr 



where the covariance matrix C is diagonal and proportional 
to the identity, /i. is a vector of the expected pixel values, 
having integrated the model over the pixel area and account- 
ing for PSF effects (which we ignore in this analysis). For this 
study, we replace the matrix summation by a continuum ap- 
proximation and integrate over the image. After some Math- 
ematica algebra, the Fisher matrix simplifies to 
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(15) 



where the rows and columns correspond to the order — 
{X,e,^Io,(^>)■ 

The correlation matrix of the parameters is formed from 
the inverse of the Fisher matrix, and is 



\ 



(16) 



From this we see that our estimate of A is completely un- 
correlated with the position angle 0, and strongly anti- 
correlated with the central surfae brightness ^o, as one might 
expect. The important result is that there is no correla- 
tion with the ellipticity e, and this is not necessarily ex- 
pected. This arises from our choice of size parameter as 
A — ^ Area/TT. The choice of the semi-major axis as the 
size parameter in Casaponsa et al. (20121 is not nearly as 



useful, as it is highly correlated with ellipticity. 

There are many effects which are not considered in this 
analysis, such as the PSF, pixelisation, a range of profiles 
and centroid errors, all of which may lead to some correla- 
tions between size and shape, but we expect on the basis 
of this calculation that these correlations will be small pro- 
vided that VArea is used for size, and we will ignore them 
in this paper. 
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Figure 2. The relative improvement in the Dark Energy FoM 
from shape plus size vs shape alone, as a function of the number 
density of galaxies, with fixed median redshift Zm = 0.9. A Euclid- 
like survey with n = 30 per square arcminute is shown by the 
vertical dashed line. 
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Figure 3. Relative improvement in marginal errors of Dark En- 
ergy parameters, when lA systematics are included approximately 
in the shape analysis. Shape alone is shown in red (outer), and 
with size added in blue (inner). The FoM is increased by over a 
factor of 4. 



4 RESULTS 

We consider a 15000 square degree survey similar to that 
proposed for the ESA Euclid mission. We assume a redshift 
distribution n{z) oc z'^ exp[—{lAlz/ ZmY'^], with a median 
redshift Zm ~ 0.9 and a mean number density fi = 30 per 



(Joachimi & Bridle 2010 Kirk et al. 20121. In the case of 



square arcminute. We assume a dispersion in InA of 0.3 ( Shen 
et al.||2003| [Ferguson et al.||2004| , and ae = 0.3 . We con- 



sider tomography with 10 bins between redshifts and 2, 
with equal numbers per bin. We compute the lensing poten- 
tial power spectrum for each bin using CAMB to compute 
the matter power spectrum, and vary the following cosmo- 
logical parameters: Q,b, ^c, ^a, h, wq, Wa, ris, 10^ A, being, 
respectively, the density parameters in baryons. Cold Dark 
Matter and Dark Energy, the Hubble parameter in units of 
100km s~^ Mpc^^, the Dark Energy equation of state pa- 
rameters (p/ p — Wo + Wa{l — a) , where a is the scale factor), 
the scalar spectral index, and the amplitude of fluctuations. 
We have not included a number of effects, such as intrin- 
sic alignments or photometric redshift errors in this proof- 
of-concept study, but the main interest here is the relative 
change in the errors when size magnification is added, rather 
than in absolute values. We show this two ways: firstly by 
showing the marginal errors of pairs of parameters, in Fig.[l] 
and secondly by computing the Figure of Merit (FoM) for 
Dark Energy, defined to be the inverse of the area/vr of the 
la contours of the expected likelihood in the wo,Wa plane, 
marginalised over all other parameters. This is shown as a 
function of n in Fig. [2] For n = 30 the FoM is increased from 
293 to 492, an improvement of 68%. 



size, there is an anticorrelation between size and luminosity, 
which reduces the size magnification, because it is accom- 
panied by fiux magnification, which brings in less luminous 
and hence smaller galaxies into the sample. This reduces the 
effect, dependent on the slope of the mean size-luminosity 
threshold. This depends on the mean size-luminosity rela- 
tion, X oc . Estimates for P vary slopes from 0.3 even 
up to unity ( jBernardi et al.||2012| ). Note that the effect is 
much diluted if the sample extends below L*, as the ad- 
ditional sources brought in by flux magnification are then 
a small proportion of the total, and the slope of the mean 
size-luminosity threshold is small. At two magnitudes below 
L* (the limit of Euclid at z ~ 1.8) the effect is a few per- 
cent only. An analogous effect to lA may exist in the form 
of size-size or size-density correlations. Studies differ in their 



conclusions with current data (Cooper et al.|2012 [Papovich 


et al. 


2012||Maltby et al."2010!IRettura et al.|2010||Cimatti 


et al. 


2012 


Park & Choi 2009), and this will need careful 



study. We will present a full study of size-shape weak lensing 
with systematics included in a later paper, but given that 
the size-magnitude effect is likely to be much smaller than 
the effects of lA, we expect the improvements presented here 
to be rather conservative. To get a rough idea, increasing CTe 
by a factor of two to approximate crudely the effect of lAs, 
and conservatively increasing crinA / by 10% for the size- 
luminosity correlation, we flnd a relative improvement in the 
FoM by a large factor of 4.2 (Fig.[3|. 



4.1 Systematic Errors 

Both size magnification and shear are subject to systematic 
errors. In the latter case, a major source is intrinsic align- 
ments (lA; e.g., [Heavens et al.|2000l [Hirata fc Seljak|2004 |. 
This can be converted into a statistical error if a flexible 
model is adopted ( [Kirk et al.[[2012| ), where marginalising 
over the lA parameters increases the Dark Energy equation 
of state errors by a factor of 2-3. Some of the lost FoM 
can be recovered with clustering information, leading to a 
degradation of Dark Energy errors by a factor of about 2 



5 CONCLUSIONS 

In this paper we have shown that adding size measurement 
to cosmic shear analyses can lead to very signiflcant im- 
provements to the Dark Energy FoM of a weak gravitational 
lensing survey. Ignoring systematics we flnd that the im- 
provement is about 68%, and we argue that much higher 
gains of even a factor 4 may be achievable when systematic 
effects are marginalised over, but this requires more detailed 
study. The full gains can be achieved if the errors from size 
and shape are uncorrelated, and we have shown that for ex- 
ponential proflles this can indeed be achieved, provided the 
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square root of the area of the source is used as the measure 
of size. We expect that for more general galaxy profiles and 
in the presence of PSF effects etc the correlation would be 
small, but non-zero. 
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